双时间步法在二维浅水方程求解中的应用          下载免费PDF全文

喻海军  黄国如  武传号 《水科学进展》2014,25(4):542-549

针对传统算法效率低的问题,将隐式双时间步法应用于求解二维浅水方程,建立了非结构网格下高效的有限体积模型。在应用双时间步法时,虚拟时间层中的定常问题采用高效的隐式LU-SGS（Lower-Upper Symmetric Gauss-Seidel）方法进行迭代求解。通过模拟计算4个典型算例以及与传统显式算法进行比较,对模型精度、效率及处理实际问题能力进行检验,分析了时间步长、内迭代次数对模型性能的影响。结果表明,双时间步法放宽了稳定性对时间步长的限制,时间步长可取到显式格式10倍以上,计算耗时减少了50%以上,模型具有良好精度与适应性,具有较好的推广应用价值。  相似文献   

关于二阶伴随模型的理论研究   总被引：9，自引：2，他引：7  

韩桂军  何柏荣  马继瑞  李冬 《海洋学报》2000,22(3):15-19

Hesse矩阵-目标函数关于控制变量的二阶偏导数形成的矩阵,在变分数据同化过程中以及敏感性分析等方面起着重要的作用;它可以通过建立数学模型的一阶和二阶伴随模型求得.以浅水方程模型为例,利用泛函的Gâteaux微分和Hilbert空间上伴随算子的概念,导出了一阶和二阶伴随模型并由此得到Hesse矩阵.改进了Zhi Wang等(1992)建立的二阶伴随模型理论.  相似文献   

Velocity predictions for shoaling and breaking waves with a Boussinesq-type model     

Franoise Ozanne  Andrew J. Chadwick  David A. Huntley  David J. Simmonds  John Lawrence 《Coastal Engineering》2000,41(4)

In this paper, the performance of a 1-D Boussinesq model is evaluated against laboratory data for its ability to predict surf zone velocity moments. Wave evolution over a plane beach and a complex bathymetry both extending into the surf-zone is examined for six cases. For the plane beach, these comprise two cases, a spilling and a plunging cnoidal wave. For the complex bathymetry, these comprise four cases of longer and short wavelengths (spilling and plunging breakers), with regular and irregular periodicity. The model evaluation places emphasis on parameters of the wave field that could be used for the prediction of sediment transport; orbital velocity, undertow, velocity skewness, kurtosis and asymmetry. It is found that, despite an overestimation of the depth-averaged horizontal velocity in the regular waves cases, the predicted higher order velocity moments and undertow are in good agreement with the laboratory data. A bispectral analysis demonstrates that the nonlinear transfers of energy amongst the low order harmonics are well reproduced, but energy exchanges with the higher harmonics are less well predicted. As a result, the model handles velocity moments better in the shorter wave tests than in the long wave cases where triad interactions are stronger. Of the four parameters describing wave breaking, the model behaviour is most sensitive to the critical wave front slope φ_B, especially with regard to velocity skewness and kurtosis predictions. It is also found that increasing the thickness of the surface roller for the case of plunging breakers improves the model's performance.  相似文献   

利用伴随法优化非线性潮汐模型的开边界条件Ⅰ．伴随方程的建立及“孪生"数值试验   总被引：1，自引：1，他引：0  

韩桂军  马继瑞  刘克修  李冬 《海洋学报》2000,22(6)

本研究采用基于最优控制理论的伴随法,把观测资料同化到陆架海域潮汐数值模型中去,优化开边界条件,提高数值预报的精度．潮汐模型的控制方程为考虑平流项、非线性底摩擦和侧向涡动粘性项的非线性浅水方程组;采用Lagrange乘子法建立了伴随模型．研究分两部分：第一部分即本文,建立非线性浅水方程模型的伴随方程、给出目标函数的梯度,并实现“孪生”数值试验;第二部分另文给出．  相似文献   

地形与流对水平无旋浅水波的影响   总被引：1，自引：0，他引：1  

孙明光  高虎山 《海洋工程》2000,18(2):62-67

本文在水平无旋及Boussinesq假设之下,导出了水面变化与水平速度场耦合方程组,以及相应的压力与垂向速度解析表示式。通过数值方法求得水面变化及某一特定水深的水平速度分布之后,压力分布及其余水深的速度分布即可由简单计算得到。由色散关系式可知,不同水深的长波色散关系在O(ε)近似之下是等同的。粘性的存在会使波高随时间的增加而衰减,但粘性与底斜率的耦合又可能使波高增长,形成不稳定;计算及分析说明,当同向流增大时,波速增大,波长增大,波高会增加;而水深减少,使波速减少,波长缩短,振幅增加。  相似文献   

不规则波Boussinesq型方程的造波,消波和反射   总被引：5，自引：1，他引：4  

朱良生  洪广文 《海洋工程》2000,18(4):43-48

对前人提出的造波、消波和反射边方法分析表明,其方法是极浅水波近似,不适用于任意水深的水域,本文就任意水深变化Ｂｏｕｓｓｉｎｅｓｑ型方程,提出了不规则波新的造波原理、方法和消波边界及部分反射边界波动方程,试验表明,本文提出的造波、消波和反射方程有效而可靠的。  相似文献   

Second kind Chebyshev operational matrix algorithm for solving differential equations of Lane–Emden type     

《New Astronomy》2013

In this paper, we present a new second kind Chebyshev (S2KC) operational matrix of derivatives. With the aid of S2KC, an algorithm is described to obtain numerical solutions of a class of linear and nonlinear Lane–Emden type singular initial value problems (IVPs). The idea of obtaining such solutions is essentially based on reducing the differential equation with its initial conditions to a system of algebraic equations. Two illustrative examples concern relevant physical problems (the Lane–Emden equations of the first and second kind) are discussed to demonstrate the validity and applicability of the suggested algorithm. Numerical results obtained are comparing favorably with the analytical known solutions.  相似文献   

A two-dimensional finite volume hydrodynamic model for coastal areas: Model development and validation     

A.H. Daoud  K.A. Rakha  A.G. Abul-azm   《Ocean Engineering》2008,35(1):150-164

A two dimensional implicit finite volume scheme for solving the shallow-water equations is developed. The effects of the Coriolis force, surface wind stress, and waves are included. A non-uniform rectilinear forward staggered grid is used with Cartesian coordinates. The time integration is performed using the Euler implicit technique. The convective flux is treated using the deferred correction method. The viscous terms are discretized using a second order central difference approximation. The SIMPLE (Semi-Implicit Method for Pressure-Linked Equations) algorithm is used for coupling the velocity components and the water elevation gradient for the water level correction. The system of equations is solved sequentially using the Strongly Implicit Procedure (SIP). To simulate wave driven current, a phase averaged wave model is used first to simulate wave transformation and calculate radiation stresses. The performance of the developed model is validated for different sources of external forces and different combinations of boundary conditions. The validation cases include tidal circulation in a harbor and wave induced currents behind a breakwater parallel to the coastline. The model is finally applied to simulate the flow pattern in a closed artificial lagoon and along the coastline near Damietta Port located along the Northern coast of Egypt. Results of the developed model agree well with the published results for the considered cases.  相似文献   

附有限制条件的间接平差秩亏时解法初探     

秦永宽  黄声享  张书毕 《海洋测绘》2009,29(1):7-9

在附有限制条件的间接平差的应用中，若所列观测方程中实际用到的未知参数的个数少于所选参数的总数u时，就会出现法方程系数阵秩亏的问题。为了解决此问题，在理论公式推导的基础上，给出了一种新的解决方法，并结合算例进行了验证。  相似文献   

Parameter estimation for basin-scale ecosystem-linked population models of large pelagic predators: Application to skipjack tuna   总被引：2，自引：1，他引：1  

Inna Senina  John Sibert  Patrick Lehodey 《Progress in Oceanography》2008,78(4):319-335

A Spatial Ecosystem and Population Dynamic Model (SEAPODYM) is used in a data assimilation study aiming to estimate model parameters that describe dynamics of Pacific skipjack tuna population on ocean-based scale. The model based on advection–diffusion–reaction equations explicitly predicts spatial dynamics of large pelagic predators, while taking into account data on several mid-trophic level components, oceanic primary productivity and physical environment. In order to improve its quantitative ability, the model was parameterized through assimilation with commercial fisheries data, and optimization was carried out using maximum likelihood estimation approach. To address the optimization task we implemented an adjoint technique to obtain an exact, analytical evaluation of the likelihood gradient. We conducted a series of computer experiments in order to (i) determine model sensitivity with respect to variable parameters and, hence, investigate their observability; (ii) estimate observable parameters and their errors; and (iii) justify the reliability of the computed solution. Parameters describing recruitment, movement, habitat preferences, natural and fishing mortality of skipjack population were analysed and estimated. Results of the study suggest that SEAPODYM with achieved parameterization scheme can help to investigate the impact of fishing under various management scenarios, and also conduct forecasts of a given species stock and spatial dynamics in a context of environmental and climate changes.  相似文献